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We report the a.c. susceptibility study of DyiTb2-a;Ti207 with x £ [0, 2]. In addition to the 
single-ion effect at Ta (single-ion effect peak temperature) corresponding to the Dy^+ spins as that 
in spin ice Dy2Ti207 and a possible spin freezing peak at Tf {Tf < 3 K), a new peak associated 
with Tb^^ is observed in Xac{T) at nonzero magnetic field with a characteristic temperature T* 
{Tf < T* < Ts). T* increases linearly with a; in a wide composition range (0 < a; < 1.5 at 5 
kOe). Both application of a magnetic field and increasing doping with Dy'^^ enhance T*. The T* 
peak is found to be thermally driven with an unusually large energy barrier as indicated from its 
frequency dependence. These effects are closely related to the crystal field levels, and the underlying 
mechanism remains to be understood. 
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I. INTRODUCTION 

Geometrically frustrated magnetic materials,—"— which 
are ordered systems with specific geometry that accom- 
modates closely competing spin-spin interactions and 
therefore frustration, attract extensive interests because 
of their novel ground states. The pyrochlores^ pro- 
vide a variety of geometrically frustrated magnetic ma- 
terials, including spin liquid,- spin ice^'^ and spin 
glasSf^"— where magnetic exchange and dipolar inter- 
actions of the nearest-neighbor spins dominate the low- 
temperature magnetic property. Among these, two py- 
rochlores are of particular interests, spin ice Dy2Ti2 07 
(DTO), HozTizOr and spin liquid TbzTiaOy (TTO), 
which have completely different ground states while shar- 
ing essential similarities. 

For the spin ice compound, spins upon each tetra- 
hedron adopt the two-in-two-out configuration in the 
ground state, which is a direct analogy to the two-short- 
two-long proton bond configuration in water icci^ This 
results in the zero point entropy. Strong crystal field 
(CF) splitting induces uniaxial anisotropy of Dy^"*" spins 
giving rise to the Ising-like ground state doublet, in DTO, 
with spins pointing along the local (111) axis, which lies 
more than 200 K below the first excited statej^^rii Novel 
spin dynamics has been observed. The single-ion effect 
at 15 K and spin freezing into the ice statei^"— (strictly, 
spins falling into compliance with the ice rule and freez- 
ing into disordered states) below 2 K are found to be 
responsible for the low-temperature dynamic properties. 
Moreover, it has been predicte d^^i^" and most recently 
verified experimentallj«2i"— that the dipole-dipole inter- 
actions give rise to magnetic monopoles, therefore a new 
wave of reconsideration of the equivalency of electricity 
and magnetism is invoked. 

The spin liquid compound TTO has Ising type spins 
similar to DTO but with a much smaller gap of 18.7 K 
between the ground and excited spin states^ii^ which 



is one order of magnitude smaller than that of DTO. 
TTO remains in spin liquid state down to 50 mK de- 
spite a short-range antiferromagnetic (AFM) orderi ^^i^^ 
With a Curies- Weiss temperature of -14 K (Ref<^), the 
absence of an ordered ground state is unusual and has 
been attributed to the delicate balance between the AFM 
exchange and ferromagnetic (FM) dipolar interactions, 
which then places TTO right at the phase boundary ac- 
cording to the phase diagram by Hertog and Gingras.— 
This delicate balance, however, is vulnerable to exter- 
nal perturbations, such as pressure and magnetic fields. 
For example, TTO shows an AFM ground state with 
spin-liquid-like fiuctuations under high pressurCf^S a long 
range order with spin wave excitations^^ and a spin-ice- 
like order with the A: = propagation vector— under high 
magnetic fields along the [110] axis. It is also proposed 
that quantum fluctuations arc responsible for the lack of 
an ordered ground state, and TTO is argued to be in a 
quantum mechanically fluctuating spin ice statei^i^ 



It is therefore interesting to ask how the ground state 
and magnetic interactions evolve if both Dy and Tb ions 
are introduced into the rare earth (RE) sites. In a brief 
reporl*^, both spin ice and liquid signatures are found 
in the ground state of Ho2;Tb2-a;Ti207. Here, using the 
a.c. susceptibility study on Dy;rTb2-a;Ti207 (DTTO), 
we observed a new peak in Xac{T) at the characteristic 
temperature T* {Tf <T* <Ts). T* is found to increase 
linearly with Dy composition x in a surprisingly wide 
range {0 < x < 1.5) for 77 = 5 kOe. The phase diagram 
of the composition dependent T* is obtained. The T* 
peak is found to be thermally driven with an unusually 
large energy barrier as indicated from its frequency de- 
pendence. The observed peaks are closely related to the 
CF, and possible origins are discussed using the compo- 
sition and the fleld dependence of a.c. susceptibilityi^ 
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FIG. 1: (color online) The lattice constant (cubic phase, a = 
b = c) of of Dya;Tb2-iTi2 07 as a function of x show a quasi- 
linear decrease. (Upper inset) X-Ray diffractions show single 
phase (pyrochlore structure) of all samples, (lower inset) Ex- 
panded XRD patterns show the peak shift with a;(Dy), the 
full width at half maximum does not show increase, demon- 
strating the absence of phase separation. 

II. EXPERIMENTAL METHODS 

Polycrystalline DTTO with x £ [0,2] were prepared 
by the standard solid-state reaction methodj^ii^. The X- 
ray diffraction data show that all samples are of single 
phase with a pyrochlore structure. The lattice constants 
(cubic phase, a = h = c) are shown in Fig. [T] For 
the two samples at the composition boundaries, lattice 
constants of 10.15(1) A for TTO {x = 2.0) and 10.12(1) 
A for DTO (x = 0) are in good agreement with pre- 
vious reports^i^ Curie- Weiss temperatures by fitting to 
the d.c. susceptibility between 10 and 20 K are consis- 
tent with values reported)^ As there are two spin species 
on the RE sites, it is important to address the issue of 
possible phase separation. As shown in the lower inset of 
Fig. [1] the full width at half maximum of the diffraction 
peaks does not show obvious increase with different a;, 
which indicates that the system is free from phase sepa- 
rations. The quasi-linear decrease of the lattice constants 
with increasing x further confirms the samples' quality. 
D.c. and a.c. susceptibilities are measured by a Quantum 
Design SQUID magnetometer and a physical properties 
measurement system (PPMS) using an excitation field 
of Hac = 10 Oe at frequency / (10 f ^ 10000 Hz) 
respectively. 



III. RESULTS AND DISCUSSIONS 

In Fig. [21 we take sample x = 1.0 as an example and 
show the real and imaginary parts of the a.c. susceptibil- 
ity at varying magnetic fields. At = 0, x' {T) increases 
monotonically with an increasing slope upon cooling for 
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FIG. 2: (color online) The real (upper panel) and imaginary 
part (lower panel) of a.c. susceptibility of Dyi.oTbi.oTi207 
for several fields. Pointed by the arrows are: Tf peak (blue 
arrow), T* peak (red arrow) and Ts peak (green arrow) de- 
noting the low temperature freezing peak, the peak associated 
with Tb"*^, and the single- ion peak associated with Dy'^"'" re- 
spectively, as discussed in Fig[3]and corresponding text. 



all excitation frequencies in study (/ = 1, 2, 5, 10 kHz), 
which signals a canonical paramagnetic behavior.— 
The imaginary part 'x"{T) in the lower panel shows a 
sharp increase as well, even when approaching our lowest 
available temperature 2 K. Therefore there should exist 
a strongly dissipative process at T < 2 K. The zero field 
a.c. susceptibility is qualitatively similar to that of TTO 
(Ref.^^) for all DTTO samples in study. The magnitude 
of x'(r) at iJ = grows with increasing a;, presumably 
due to the larger moment of Dy'^"'" spins. As the field 
increases to 5 kOe, a clear dip in x'{T) appears around 
8.5 K, leading to a local maximum in x'(r) as well as a 
correlated peak in x"(2^) dictated by the Kramers-Kronig 
relation. It is also notable that there is an emerging max- 
imum on the right shoulder of the peak in x"(T), which 
becomes prominent as the field increases. At H = 10 
kOe, x'(^) shows two clear dips at 12 K and 17.5 K, 
thereby two local maxima appear in x"{T). As the fre- 
quency increases, the separation between the two local 
maxima increases (Figl2]c"). 

The two spin-freezing-like peaks cannot be attributed 
to a collective spin glass transition, although the random 
occupancy of rare earth sites by Dy^^"*" and Tb^"*" is a 
likely source of spin glassiness. In contrast to the spin 
glass transition, the d.c. susceptibility of DTTO shows 
simply a Curie- Weiss behavior, more strictly, a superpo- 
sition of Curie- Weiss and Van Vleck termsj^ii^s^^ with no 
freezing transition. Secondly, the application of a mag- 
netic field enhances the characteristic temperature of the 
two freezing-like peaks, as opposed to the case of spin 
glass where the transition temperature is suppressed.— 
Furthermore, the interactions for Dy-Dy and Tb-Tb on 
the ground state doublet in DTO and TTO cases are only 
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FIG. 3: (color online) The real (upper panel) and imaginary part (lower panel) of a.c. susceptibility of Dy2,Tb2-a;Ti207 {x = 
0, 0,5, 1.25, 1.5 and 1.8) at // = 5 kOe. Pointed by the arrows are: low temperature freezing peak (blue arrow), the single-ion 
effect peak associated with Dy^^ (green arrow) and the T* peak (red arrow). Insets in the a; = plots are the real (upper 
inset) and imaginary part (lower inset) of a.c. susceptibility of Dya;Tb2-i:Ti2 07 with x = at H — 10 kOe. 



of the order of 1 K, much lower than the temperature of 
the observed peaks. Therefore the freezing-like peaks are 
most likely not due to spin glass transitions. With higher 
magnetic fields {H = 10 kOe data shown in Fig. IH^b) 
and higher field data not shown here), the thermal spin 
contribution is suppressed at low and high temperatures, 
leaving only a broad peak in x'{T) as observed in TTO 
and DTOi^ Other samples with 0.5 < cc < 1.5 show 
qualitatively similar behaviors. The appearance of the 
two freezing-like peaks as well as their strong frequency 
dependence are reminiscent of the single-ion effect peak 
in pure DTO, i.e., x = 2.0, and will be discussed in detail 
by examining the composition dependence below. 

Fig. |3] shows the a.c. susceptibility of five representa- 
tive samples with a; = 0, 0.5, 1.0, 1.5 and 1.8 at _ff = 5 
kOe. For x = 0, there appears to be a peak below 3 K 
in x'iT). The imaginary part of the sample with x — 
is clearly different from other doped samples regarding 
their frequency dependence. This is presumably related 
to the fact that the exchange coupling and the energy gap 
between narrowly spaced CF levels of TTO are of compa- 
rable energy scale and therefore both have nonnegligible 
contributions to the susceptibilityi^i^ For samples with 
< cc < 0.5 (data not shown here), the crystal field con- 
tribution increases progressively as x increases, manifest- 
ing as the transition to the behavior of the sample with 
x = 0.5. The sample with x = 0.5 shows a single peak in 
x"(T) (T* indicated by the red arrow). Additional peaks 
(Tf and Tg indicated by the blue and green arrows) start 
to emerge for x — 1.0, and become more pronounced in 
X — 1.5 and 1.8 samples. These additional peaks for 
X ^ 1.0 have the same origin as the two peaks in pure 
DTO. In pure DTO, two peaks are observed in x"{T), 
a high temperature peak corresponding to the single-ion 
effect, and a low temperature peak associated with the 
spin freezing dictated by the ice-rule ground stateJ^"— 



For the sample with x ~ 1.8, the position of the high tem- 
perature peak (indicated by the green arrow) are in good 
agreement with the single-ion peak in DTO. Moreover, 
the frequency dependence of the high temperature peak 
and the extracted energy scale corresponding to the crys- 
tal field splitting are qualitatively similar to that of pure 
DTO, as discussed below. We therefore identify the high 
temperature peak as the single-ion peak of Dy'^"'" spins. 
The low temperature peak (indicated by the blue arrow) 
also possesses a similar characteristic temperature with 
the spin freezing peak in DTOfi^ the corresponding peak 
in x"iT) should appear below 2 K, the limit of PPMS, 
thereby only a rapid increase in x"(T) upon cooling be- 
low 5 K is observable, and a direct comparison with the 
DTO data is not allowed. Following the convention of 
previous reports, we define the characteristic tempera- 
ture Ts as the high temperature peak position (green ar- 
row) in x"{T) and similarly Tf for the low temperature 
peak (blue arrow). 

The most important new feature in all samples stud- 
ied, is the presence of a third peak (indicated by the red 
arrow) between Tf and Ts. We define T* {Tf < T* < Ts) 
as the characteristic temperature for this third peak. T* 
peak should root in the Tb'^"'" spins as it exists at regions 
with small x (including x = at iJ = 10 kOe, see the 
insets in Fig. [3]) , where the peaks associated with Dy'^^ 
disappear and the effect of Dy^+ should play a minor role. 
Combined with the susceptibility of Dyi.oTbi.oTi207 at 
different fields shown in Fig. [2l it is interesting to note 
that the evolution of the a.c. susceptibility with increas- 
ing X (Fig. 131) is similar to that with constant x but 
increasing magnetic field (Fig. [5]). In both cases, the in- 
crease of X or magnetic field leads to the appearance of 
the T* peak and subsequently the Ts peak. 

In order to investigate the origin of the T* peak, we 
studied a series of samples with x G [0, 2], and track T* 
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FIG. 4: (color online) Contour plot of x"{T) of samples with different x ai H = 5 and fO kOe with / = 10 klfz. Pink circles 
(blue squares) denote T* (Ts). Dotted lines are guide to the eye where the peak positions are not well resolved. 



and Ts in all samples when observable. A global view 
of x"{T) with different x is shown in the contour plot in 
Fig. 01 It is clearly seen that the T* peak exists in a wide 
range of x. Most interestingly, T*{x) increases linearly 
with increasing x in a wide range of x, both at 5 kOe and 
10 kOe. At _ff = 10 kOe, the low temperature upturn is 
strongly suppressed compared with the 5 kOe case, as 
mentioned above. The linearity of T*{x) persists in a 
wide x range, from x = to 1.5, obvious deviation from 
the linearity occurs when x > 1.5. The characteristic 
peak temperature of Tg is also shown in Fig. |4l with a 
slight decrease as x becomes smaller. 

A clear deviation from linearity in T* for x ^ 1.5 is 
also an indication for its Tb'^"'' origin. The pyrochlore 
structure of DTTO consists corner-shared tetrahedrons, 
by assuming that Dy^+ and Tb^+ spins in DTTO oc- 
cupy the lattice randomly therefore uniformly, there are 
two critical values of a;(Dy) where large change in local 
environment (i.e. tetrahedron that consists of only one 
spin species: all Dy'^^ or all Tb'^"*") emerges: for x > 1.5 
{x ^ 0.5), tetrahedrons of Dy'^+-only (Tb^+-only) spins 
appears, while Tb'^+(Dy'^+) still lies in the uniform neigh- 
borhood coordinated with six Dy'^^ (Tb'^^) spins. The 
deviation from the linearity in T*{x) can be understood 
as a direct consequence of the critical value of x = 1.5. 
It is important to note that there is no such deviation at 
the Tb end. The presence of the deviation at a; = 1.5 as 
well as the lack of that at a; = 0.5 is a strong indication 
of the asymmetric roles of Tb'^"'" and Dy"^"*" spins in the 
origin of T* peak, consistent with our argument that the 
T* peak is associated with Tb^+. 

To see the nature of T*, we show in Fig. [Sja) the 
phase diagrams of the Ts{x) and T*{x) with different a.c. 
excitation frequencies (1, 2, 5 and 10 kHz) at iJ = 10 



kOe. The quasilinear increase of T*(x) is preserved at 
H — 10 kOe at different frequencies. The frequency de- 
pendence of both T* and Tg behaves according to the 
Arrhenius relation, as shown in Fig. [5jb). One note that 
to achieve a reliable fitting with the exponential function 
/ = foexp{—EA/kBT), data with frequencies spanning 
over several decades are necessary. In our case, how- 
ever, at frequencies lower than 1 kHz, T* and Tg peaks 
are indistinguishable, making the data at lower frequency 
range inaccessible. Thus we only focus on the data above 
1 kHz. From the fitting of / = foexp{—EA/kBT), where 
/o is taken to be the physically reasonable values of 1 
GHzfiiii^ the energy barrier Ea can be extracted. The 
extracted energy barrier Ea is shown in Fig. [Sfc), includ- 
ing Ea by fitting with /q = 0.5 and 5 GHz to verify the 
reliability of the fittings. The validity of the fitting by Ar- 
rhenius law indicates that in the presence of a magnetic 
field the spin relaxation is thermally drivenJ^ In pure 
DTO, Ea has been found to be 200 K (Ref-i^). This en- 
ergy barrier has been interpreted as the gap between the 
lowest lying crystal field levels of Dy"^+. In DTTO, it can 
be seen that Ea associated with Tg changes only slighly 
from 193 K for x = 1.25 to 200 K for x = 1.8, a value 
very close to the energy barrier in pure DTO. This again 
proves our understanding that Tg peak is associated with 
the single-ion effect of Dy'^^. On the other hand, Ea as- 
sociated with T* (attributed to Tb^"*" above) increases 
rapidly from 85 K for x = 0.25 to 154 K for x = 1.0. 
These values are unusually high, since the crystal field 
gap of ground state doublets for TTO is only 0(20 K). 
Similar fittings with data at 5 kOe cannot be performed 
since the T* and Tg peaks are not well resolved for low 
frequencies. 

Our observations including the field, composition and 
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FIG. 5: (color online) (a) At = 10 kOe, T* (empty sym- 
bols) and Ts (solid symbols) as a function of 2;(Dy) at different 
frequencies. Dash lines are guide to the eye to show the quasi- 
linear increase of T*. (b) Fittings by / = foexp{—EA/kBT) 
with fo taken to be the physically reasonable value of 1 GHz. 
(c) The extracted energy barrier Ea for T* (open symbols) 
and Ts (solid symbols) with fitting parameter fo — 0.5, 1 and 
5 GHz. 



frequency dependence of T* raise the question of its un- 
derlying mechanism. Knowing the exchange and dipolar 
interaction in both DTO and TTO (for TTO, nearest 
neighbor exchange J„„ = -0.88 K, nearest dipolar in- 
teraction D„„ = 0.8 K (Ref,^), for DTO J„„ = -1.2 
K, D„„ = 2.35 K (Ref.-'')), it is hard to imagine that 
certain combinations of interactions between Dy and Tb 
spins (i.e. Dy-Dy, Tb-Tb and Dy-Tb) can give rise to 
a characteristic temperature T* with its energy barrier 
two orders of magnitude higher than the energy scale set 
by the interactions. Instead, the large energy gap ex- 
tracted from the frequency dependence of T* suggests 
that single-ion effects may play an important role. The 
CF levels for pure TTO and DTO have been studied in 
detail previously/i^^!^^ For the DTTO samples, there 



exist two sets of CF schemes associated with Dy^+ and 
Tb'^"'" respectively, with their CF level structures deter- 
mined by the local chemical environment. The zero-field 
Ea in pure DTO is around 200 K. With = 10 kOe, Ea 
associated with Dy'^^ in DTTO also stays around 200 K 
as is shown in Fig. [DJc). Thus one may expect that the 
single-ion peak temperature corresponding to the CF 
of Dy'^^ in DTTO should be robust against the change 
in magnetic field. However, according to the a.c. sus- 
ceptibility data at zero field, e.g., Fig. [SJa' ) and (a" ), 
none of the sample with x ^ 1.5 show signs of a single- 
ion peak above 2 K, implying narrowly spaced CF levels 
similar to that in TTO. Moreover, the field dependence 
of T* and in Fig. [5] indicate a dominant impact on 
their respective CF gap by the field. The possible field 
dependence of the CF levels is novel and deserves fur- 
ther investigation. Here with the current limited data, 
we discuss possible mechanisms that may give rise to 
the field dependence. It is known that with the appli- 
cation of a magnetic field, the Zeeman splitting of the 
spin states may cause changes in the CF scheme. Take 
Dy-^^ ions as an example^^ gjJ ^ O{10)^b / Dy. With 
the magnetic field of 1 T, the Zeeman energy is gjJH ~ 
0(10) K, which is one order of magnitude smaller than 
the experimental Ea value. We therefore reason that 
there should be some interactions beyond the CF inter- 
action that contribute to the magnetic field dependence. 
Previous studies on TTO and similar geometrically frus- 
trated systems revealed the softening of the crystal field 
gap due to the presence of the spin-spin correlation, 
The spin-spin correlation can significantly alter the CF 
in TTO since its energy gap is only 0(10) K. It is there- 
fore possible that the interplay between CF interaction 
and spin-spin correlation dominate the magnetic field de- 
pendence of the single-ion peaks. At zero field, the gap 
of CF is softened due to the spin-spin correlation be- 
tween Dy3+ and Th^+ , causing the disappearance of the 
Ts and T* peak. With the presence of a field, the spin- 
spin correlation is dramatically altered, the gap of CF 
is again activated and thus the peaks appear. However, 
the spin-spin interactions are orders of magnitude smaller 
than the gap observed here, thus may not cause sizable 
softening of the gap. Another possibility is the pres- 
ence of magneto-elastic coupling^Si^i^. A giant mag- 
netostriction has been reported in TTO (Ref 4^), indi- 
cating strong coupling between spin and lattice degree of 
freedom. In Dya;Tb2-a;Ti207, it is possible that similar 
magneto-elastic coupling act in a way that the magnetic 
field may change the chemical environment of the RE 
ion, and tune the crystal field effectively. However, such 
a strong effect on crystal field (change by several tens of 
Kelvin) by a moderate magnetic field (1 T) is still sur- 
prising. The detailed mechanism remains to be explored. 

At = 10 kOe, with Ea associated with Dy^"*" at- 
taining its zero-field value of 200 K, the composition de- 
pendence of the T* peak can be understood as a conse- 
quence of the change of the crystal field with Dy'^+ dop- 
ing. The crystal field scheme is intimately related to the 
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structural parameters of the systemf24 including the lat- 
tice constants and the positional parameters of the eight 
oxygen surrounding the RE ion. Dy^+ doping leads to a 
systematic change in the structure and the corresponding 
changes in the crystal field levels. A very recent specific 
heat measurement on the same system^* shows that for 
small concentration of Dy, a high temperature peak cor- 
responding to the first excited CF level of Tb^+ moves to 
higher temperature, indicating a modification of the CF 
of Tb^^ site, which is consistent with our understand- 
ing. Also, the linearity of T*{x) can be understood as a 
consequence of the monotonic increase of Ea{x) by the 
Arrehnius relation with suitable parameters. Another 
possibility for the composition dependence of T* peak is 
the crystal field-phonon coupling as suggested in TTO,'*'* 
as Dy'^'*' doping alters the phonon spectrum and thus the 
CF level scheme. 

The aforementioned similarity between the evolution 
of the a.c. susceptibility with increasing Dy doping x 
and that with constant x but increasing magnetic field 
now has a plausible interpretation. In both cases, the 
change in crystal field levels dominates the behavior of 
the a.c. susceptibility: in the former case, the Dy doping 
changes the structure and alters the crystal field directly; 
in the latter case, the magnetic field affect the crystal field 
possibly through the magneto-elastic coupling. 

IV. CONCLUSION 

To conclude, we have systematically studied the a.c. 
susceptibility of Dya;Tb2-3;Ti207. A new peak associ- 



ated with single-ion Tb^"*" (T* peak) is observed, to- 
gether with another peak associated with that of Dy'^+ 
(Ts peak). The phase diagram of both T*{x) and Ts{x) 
is presented showing the characteristic temperature in- 
creasing linearly with x. Both Ts and T* peaks show 
strong frequency dependence, suggesting that the transi- 
tion is thermally activated. We interpret the origin of the 
energy barrier to be associated with the gap in the low- 
est lying CF levels of Dy'^^ and Tb'^"'". However, the ex- 
tracted energy barrier for the T* peak is an order of mag- 
nitude higher than that of TTO and comparable to that 
of DTO. Furthermore, despite the large energy barrier, 
0(100 K), these characteristic temperatures are strongly 
influenced by applied magnetic field. Further neutron ex- 
periments and theoretical calculations may help to fully 
understand the system. 
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